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Understanding the dynamics around rotating black holes is imperative to the suc- 
cess of the future gravitational wave observatories. Although integrable in principle, 
test particle orbits in the Kerr spacetime can also be elaborate, and while they have 
been studied extensively, classifying their general properties has been a challenge. 
This is the first in a series of papers that adopts a dynamical systems approach to 
the study of Kerr orbits, beginning with equatorial orbits. We define a taxonomy 
of orbits that hinges on a correspondence between periodic orbits and rational num- 
bers. The taxonomy defines the entire dynamics, including aperiodic motion, since 
every orbit is in or near the periodic set. A remarkable implication of this periodic 
orbit taxonomy is that the simple precessing ellipse familiar from planetary orbits 
is not allowed in the strong-field regime. Instead, eccentric orbits trace out preces- 
sions of multi-leaf clovers in the final stages of inspiral. Furthermore, for any black 
hole, there is some point in the strong-field regime past which zoom-whirl behavior 
becomes unavoidable. Finally, we sketch the potential application of the taxonomy 
to problems of astrophysical interest, in particular its utility for computationally 
intensive gravitational wave calculations. 
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I. INTRODUCTION 

When Einstein realized his General Theory of Relativity correctly yielded the preces- 
sion of the perihelion of Mercury, he had heart palpitations, later writing to Ehrenfest [jj , 
"For some days I was beyond myself with excitement." The anomalous precession of the 
perihelion of Mercury was the only astronomical observation in conflict with Newtonian 
gravity. Even taking account of perturbations from other planets, astronomers still saw 
Mercury's perihelion overshoot its Keplerian target by an extra 43" /century, exactly the 
amount predicted by General Relativity. 

Just shy of a century later, we are on the brink of another important test of General 
Relativity: the direct detection of gravitational waves. Black hole binaries may be the 
most viable candidates for a first direct detection. Consequently, and interestingly, the first 
gravitational waves observed will likely also probe relativistic orbits, which display exotic 
behaviors beyond even Mercury's simple precessing ellipse. Relativistic orbital trajectories 
need not lie in a plane. They can also exhibit so-called "zoom-whirl" behavior - an extreme 
form of perihelionprecession - whirling around the central black hole before zooming out 
quasi-elliptically 2|, |3[]. Nevertheless, even though orbits of a non-spinning test particle 
around a spinning black hole have been studied since Carter proved their integrability in 
the 1970's 4J, a language for making simple, general claims about their properties has been 
elusive. 

In this paper, we introduce a powerful taxonomy that defines the full range of orbital 
dynamics in the equatorial plane of a Kerr black hole. We use this scheme to illustrate 
behaviors that run counter even to our relativistic intuitions. In the strong-field regime, 
precessing elliptical orbits such as Mercury's are excluded. Instead, at close separations, 
eccentric orbits trace out precessions of patterns best described as multi-leaf clovers. We 
can also demarcate a region in orbital parameter space where the aforementioned zoom-whirl 
behavior is not merely prevalent but unavoidable, the size of that region increasing as the 
spin of the central black hole increases. 

Not all orbits are created equal 

Our taxonomy emphasizes a dynamically special set of orbits - the periodic orbits that 
return exactly to their initial conditions after a finite time. First widely touted by Poincare 
jfj, who suggested that the general behavior of any classical system could be gleaned from 
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a study of repeating motions, periodic orbits have played a crucial role in the treatment of 
some difficult problems in celestial mechanics, including the motions of planetary satellites, 
the long term stability of the solar system, and motion in galactic potentials. In contrast, 
periodic orbits in relativistic astrophysical systems like compact object binaries have gone 
largely unexamined 1 , typically with the disclaimer that because they have measure zero in 
the space of all possible orbits, these closed orbits merit correspondingly little attention. 
Particularly in the case of test particle motion in the Schwarzschild and Kerr spacetimes - 
the only analytically soluble relativistic orbital systems - mapping out the properties of a 
particular measure zero set has seemed unnecessary. 

By contrast, we will side with Poincare. To paraphrase Orwell, while all measure zero 
sets are equal, some are more equal than others. For instance, circular orbits receive special 
attention, even though they are also a measure zero set, because they have two special 
dynamical features: 

1. circular orbits are easy to handle, and 

2. some orbits look like small perturbations to circular ones. 

As a result, an analysis of circular orbits reveals fairly detailed information about nearby 
low-eccentricity orbits with relatively low overhead. Most relativity texts derive the famed 
precession of the perihelion of Mercury in precisely this way [l4j. Of course, since most 
orbits are not close to circular ones, circular orbits alone do not encode the entirety of black 
hole orbital dynamics. 

Remarkably, periodic orbits do. In fact they have even greater dynamical power than 
circular orbits because 

1. periodic orbits are also easy to handle, and 

2. all generic orbits look like small perturbations to periodic ones. 

The latter fact, first noted by Poincare, stems from a beautiful correspondence between 
periodic orbits and the rational numbers. The density of rationals on the number line thus 
implies a corresponding density of periodic orbits in the space of all possible orbits, so that 



While some authors 0, S, @, 0, [HI ; 0] have explored inherently relativistic periodic motion, the 
systems studied have been more mathematically informative than astrophysically descriptive. 
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FIG. 1: This figure is a preview of figure [T2l Rows 1 and 3 show exactly periodic orbits. Rows 2 
and 4 show nearby aperiodic orbits. 
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any generic orbit can be viewed as an arbitarily small deviation from some exactly periodic 
counterpart. 

The result is a highly geometric skeleton of periodic orbits in terms of which the properties 
of even generic, aperiodic orbits can be described. Figured] offers a preview of the anatomy of 
this skeleton. Ignoring for the moment the details of these orbits, all of which are explained 
in the body of the paper, we can still see at a glance the two special dynamical properties 
discussed above. Rows 1 and 3 show a set of exactly periodic orbits. Besides being visually 
elegant, those correspond to a rational number according to a scheme detailed in ^JTTJ Just 
below each periodic orbit is a generic aperiodic orbit. Notice that each aperiodic orbit 
looks like a slow precession not of an ellipse but of the periodic orbit immediately above. 
Additionally, each such orbit can be assigned an approximating rational. What's more, we 
can describe all black hole orbital dynamics by such a periodic table. 

This paper is the first in a series that realizes Poincare's dictum for equatorial orbits in 
the Kerr spacetime 2 . After defining and filling out this periodic skeleton for equatorial orbits 
(^ITl). we extract its dynamical consequences ( §1111) . and explore its potential applications to 
a host of astrophysical problems ( §IVI) . In particular, we outline how periodic orbits might 
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Although important and interesting in its own right, for simplicity of presentation we 
relegate to an appendix the Hamiltonian formalism we use to identify periodic orbits in 
the Kerr system. For the body of the paper, the reader need only know that a is the spin 
parameter for each central black hole and that every orbit is specified by the energy E and 
angular momentum L of a test-particle as measured by an observer at infinity. As detailed 
in appendix we work in units in which a, E, L, and all coordinates and frequencies are 
dimensionless. 

II. TAXONOMY OF PERIODIC ORBITS 

The goal of this section is to detail a system for indexing all closed orbits around a 
black hole with a triplet of integers (z,w,v). The reason behind the choice of symbols 



2 For clarity of exposition, we relegate the dynamically more involved nonequatorial Kerr orbits to another 
work. 
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FIG. 2: Left: Half of the (z = 2, w = 0, v = 1) periodic orbit. Right: The full (z = 2, w = 0, v = 1) 
closed orbit. The orbit has a = 0, L = 3.980393, and E = 0.973101. 

will become clear shortly. The scheme is topological, and our approach is to establish the 
connection between a given periodic orbit and its (z, w, v) label visually. We then establish 
the relationship between a given periodic orbit and a specific rational number in two ways: 
first based on topological features of the orbit, and then based on frequencies associated 
with its radial and azimuthal motions. 

Finally, there is the matter of how we know that any particular periodic orbit we refer- 
ence even exists. Answering this rather important question is the central objective of the 
dynamical section of the paper, §1111 

A. The essential taxonomy 

We have a simple topological method for identifying all of the closed orbits for a given 
angular momentum L around a given black hole (fixed by the spin a). Each periodic orbit 
traces out a finite number of leaves before closing. We will call the number of leaves z for 
zoom. 3 

Figure [2] shows a z = 2 orbit in the right panel, with a single radial cycle of the same 
orbit (from one apastron to periastron to the next apastron) shown in the left panel. Notice 
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3 Consider again the periodic orbits in rows 1 and 3 of figure [TJ An immediately recognizable integer 
associated with each such orbit is the number of leaves it traces out before returning to its starting point. 
Each such leaf corresponds to the quasi-elliptical "zooming" behavior mentioned in the introduction. 
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FIG. 3: Left: Half of the (z = 2, w = 1, v = 1) periodic orbit. Right: The full (z = 2, w = 1, v = 1) 
closed orbit. The orbital parameters are a = 0, L = 3.718679, and E = 0.966555. 

from the figure that the accumulated angle from apastron to apastron, Aip r , is 3ir. In a 
complete orbit, the accumulated angle is Aip = zA(p r = 6ir. 

This is not, however, the only kind of z = 2 orbit. Figure [3] shows another 2-leaf orbit that 
whirls around the center an additional 2n before zooming out to apastron again. In fact, 
there is a 2-leaf orbit that whirls 4-7T longer around the center before returing to apastron, 
and in general one that whirls 2ttw longer around the center before returning to apastron, 
for any positive integer w. Therefore, we will distinguish orbits by their number of whirls, 
w, as well as by z. 

Still, z and w alone are not sufficient to specify the geometric features of a periodic 
orbit. To see this, note that the successive apastra of a periodic orbit with z > 2 form 
the vertices of a regular polygon. We will label the vertices of these polygons with a third 
integer v, counting the starting apastron of the orbit as v = and increasing in the same 
rotational sense as the orbit (counterclockwise for prograde orbits, clockwise for retrograde 
orbits), as shown in figure HI Now, given any z > 2, an orbit might move from the starting 
apastron immediately to the next vertex in the polygon. Such an orbit will be labelled v — 1. 
However, an orbit with the same (z, w) might skip the next neighbor vertex. We will assign 
that orbit a v of 2. In general, a periodic orbit with a given z > 2 can skip any number 
of vertices less than z when moving between successive apastra. All orbits will therefore 
be specified by (z, w, v) where v indicates the first vertex hit by the orbit after v — 0, and 
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FIG. 4: 4-leafed orbits with no whirls. Left: Leaves are traced out in sequential order for the 
[z = 4, w = 1, v = 1) closed orbit. Right: Leaves are traced out of order for the (z = 4, w = 1, v = 3) 
closed orbit. The orbital parameters are a = 0, L = 3.834058 for both. The energy of the leftmost 
orbit is E = 0.979032 and the energy of the rightmost orbit is E = 0.979842. 

where v has the range 

l<v<z-l (1) 

The orbit on the left of figure H] for instance is a (4, 1, 1) while the orbit drawn on the right 
is a (4, 1,3). We can still use eqn. ([T]) for the z = 2 leaf orbits, despite the fact that two 
points do not trace out a polygon, so that v = 1. Following this rubric, the orbit of figure [2] 
is a (2, 0, 1) and that of figure [3] is a (2, 1, 1). 

Single leaf orbits, such as those shown in figure El need separate discussion. Instinctively, 
we want to assign them a z value of 1, but then the v restricted by the range in ([1]) would 
be undefined. We can handle this problem by assigning such orbits v — 0, since successive 
apastra for these orbits are actually the same single apastron. That would lead us to modify 
the allowed range of v to 

1 < v < z — 1 , if z > 1 

v = o , if z= 1 . (2) 

There is another degeneracy to address. For a given w, some (z, v) pairs describe the 
same orbit. For instance, the (4, 1, 2) orbit closes after only two leaves and is identical with 
the (2, 1, 1) with the same orbital parameters. We remove this degeneracy by requiring that 
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FIG. 5: Left: The full [z = l,w = l,v = 0) closed orbit. The orbital parameters are a = 0, L = 
2.714326, and E = 0.932703. Right: The full [z = 1, w = 2, v = 0) closed orbit. The orbital 
parameters are a = 0, L = 3.535534, and E = 0.948491. 

z and v be relatively prime. The allowed v are then 

1 < v < z — 1 , if z > 1, z mod v = 

v = , if z = 1 . (3) 

Having pruned those orbits out of any counting, the (z, w, v) label now uniquely specifies 
the topological features of a closed orbit. 



B. Periodic orbits and rational numbers 



Our (z,w,v) taxononmy naturally forges the relationship between rational numbers and 
periodic orbits. Note first that any positive rational number q can be written in the form 

q = s + - , 4 

n 

where s > is the integer part and m/n is the fractional part, with m and n relatively 
prime integers satisfying 

1 < m < n - 1 (5) 

Of course, these are exactly the conditions that w, v and z satisfy. Since every periodic orbit 
corresponds to a (z, w, v) set in our scheme, we can therefore associate a rational number 

q = w + - (6) 

z 
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to every periodic orbit. This association reflects the physical observation that for a periodic 
orbit, the accumulated azimuth between successive apastra must be given by 

Aip r = 2tt ( 1 + w + -) = , (7) 



where the total accumulated angle Aip in one full orbital period is zAtp r . 

There is another sense in which we can associate a rational number to a periodic orbit. 
Every eccentric equatorial orbit has 2 associated orbital frequencies: a radial frequency 

2tt 

= TfT , (8) 
1 r 

where T r is the (coordinate) time elapsed during one radial cycle (not the total period of an 
orbit), and an angular frequency 

* r d Ji dt = %l (9) 

J-r Jo dt l r 

corresponding to the time-averaged value of dip/dt over one radial period. For a generic 
orbit, the ratio of those frequencies can be arbitrary, but for a periodic orbit, it must satisfy 

^ = ^ = i + tu + H = i + g . (10) 

The periodic orbits, then, are those whose fundamental orbital frequencies are rationally 
related. 

It is worth noting that the orbital frequencies u r and are the same frequencies that arise 
in an action-angle Hamiltonian formulation of the dynamics 23j. We draw the connection 
explicitly in appendix IB1 



C. Circular orbits 



Circular orbits, while strictly periodic, do not fit into our indexing scheme. This is not 
surprising, since circular orbits close compulsorily rather than by the tuning of any orbital 
parameters. Put in other terms, a circular orbit has only a single rotational frequency 
Up] there is no additional librational frequency to which it must relate rationally in order to 
close. Circular orbits exhibit a periodicity of a different flavor, and no association to rational 
numbers is required. They are, in a sense, periodic by default. 

Perhaps unexpectedly, there is nonetheless a natural way to fit the stable circular orbits 
into our taxonomy. Consider, for a given black hole spin a, all the orbits with a given angular 
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momentum L. For such a fixed a, L pair, the frequencies and u r vary continuously 4 with 
eccentricity. Not surprisingly, the zero eccentricity limit of is just the dip/dt of the stable 
circular orbit. Additionally, as detailed in appendix [Cj the zero eccentricity limit of u r is 
not zero but rather the frequency of radial oscillations for small perturbations of the circular 
orbit. In this limiting sense, then, we can assign every stable circular orbit an effective 



and thus an effective (z,w,v) triplet. 

As shown in Appendix [Cj q c grows monotonically with decreasing r c , the radius of the 
stable circular orbit. In fact, it turns out that for a given a, every positive rational q 
corresponds to some stable circular orbit, with q c — > as the circular radius r c — > oo and 
q c — > oo as r c — > risco, the radius of the innermost stable circular orbit. At low r, q c will 
eventually go above 1 so that w c > 1. Additionally, the corresponding z c will sometimes 
be small (z < 10). Both these facts will leave their footprint on orbits in the strong-field 
regime. 



Our taxonomy focuses on the set of measure zero closed orbits. A generic orbit will 
experience an accumulated angle that is irrational so the orbit precesses and never closes. 
However, since the rationals are dense on the number line, any generic orbit can be ap- 
proximated arbitrarily closely by some rational. Therefore, although it seems that we are 
describing a special set of extremely rare orbits, we might as well be describing every orbit. 
The potential power of this taxonomy lies in this simple observation. 

To illustrate the approximation of a generic orbit by periodic ones, consider the (3, 1, 1) 
closed multi-leaf orbit on the left of figure El There are aperiodic orbits which lie close to 
this orbit. For instance consider an orbit that accumulates a slightly greater angle in one 
cycle from apastron to apastron, w + v/z = 1 + 1/3 + 5 where the irrational 5 << 1. It's 
three leaf pattern will precess, never closing. The precessing orbit could be approximated 
by the periodic orbit (3, 1, 1). The smaller the precession, the better the approximation. Or, 




(11) 



D. Precession of periodic orbits 



4 In fact, as we'll see in Sections MI Al and IIII B[ they vary monotonically with energy and eccentricity at 
fixed a, L. 
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FIG. 6: Left: The closed (z = 3,w = l,v = 1) orbit, right: The closed [z = 300, w = l,v = 103) 
orbit looks like a precession of the (z = 3, w = 1, v = 1) orbit. Only 6 of the 300 leaves are shown. 
The first segment from apastron to apastron is emphasized in bold. The orbital parameters are 
a = 0, L = 3.834058 for both. The leftmost orbit has energy E = 0.979304 and the rightmost orbit 
has energy E = 0.979331. 

we could do better. For the sake of argument, let's approximate the irrational drift by the 
rational 5 « 1/100 and write w + v/z = 1 + 1/3+1/100 as w + v/z = 1 + 103/300. We can 
now approximate the precessing orbit by the (300, 1, 103) multi-leafed orbit on the right of 
figure EJ Our rational approximate is a 300-leaf orbit. Since v = 103, it skips 102 leaves in 
the pattern each time it moves out to apastron so that it appears to trace out a precessing 
three leaf pattern. The periodic approximate does actually close but only after tracing all 
300 leaves. 

We emphasize that the approximation of generic orbits by a nearby periodic is a reality 
foisted on us by finite precision, both numerically and observationally. Every computer 
program truncates numbers to some finite precision and in so doing explicitly outputs a 
rational. Every measurement performed delivers an observable to some finite precision and 
in so doing explicitly delivers all data as rationals. Likewise, every calculated aperiodic 
black hole orbit will necessarily be indistinguishable from some periodic one. We have 
no choice about this. The taxonomy embraces that fact and allows us to estimate the 
rational approximate deliberately rather than inadvertently through the restrictions of finite 
precision. 
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E. Not all orbits allowed 



An obvious question to ask is whether all periodic orbits are allowed in every Kerr system. 
In other words, is every rational in the set < q < oo allowed? The answer will turn out to 
be a very interesting no. 



III. BLACK HOLE ORBITAL DYNAMICS 



Indeed, not all periodic orbits are allowed in a given Kerr system. We now show that the 
rational numbers occur in the range 

Qc < q < Q max i 

(12) 

bounded below by the q c of the stable circular orbit and above by the q max of the maximum 
energy bound orbit. 5 Periodic orbits corresponding to all rational numbers in this range 
populate the phase space. The limits vary not only for different black hole spins, but also 
for different values of the constant angular momentum L, as we will demonstrate in §111 Al 
and gnUl 

Equation (fl2"|) turns out to be profound. The crucial insight that allows us to determine 
the admissable periodic orbits is the observation that A<p r , and therefore the associated 
rational q = w + vj z, increase monotonically with energy. We extract the dynamical conse- 
quences of this fact, contrasting the weak and strong- field behaviors. 



A. The Schwarzschild zoo 



To appreciate the significance of observation ffl2l) . we first consider Schwarzschild (a = 0) 
orbits because they admit a formal effective potential formulation that lends clarity to the 
discussion. Schwarzschild orbits can be described as one-dimensional motion in an effective 
potential Q: 

\r 2 + Kff = e eS (13) 



where 

-E V ---- — -- 

2 ' eff 2 r 2r 2 r 



s cS = -E, VcS = --- + — 2 --^ . (14) 



5 To be definite, we mean non-plunging, bound orbit. 



11 



■J- 




le+08 



le+09 



le+10 



FIG. 7: The effective radial potential V e Q (r) for very large r (plotted on a In r scale) in the vicinity 
of solar system values. The horizontal line at V q q = 1/2 corresponds to particle energy E = 1. 

Since the shape of the potential is fixed by the value of L, a given V e s snapshot is tantamount 
to a snapshot of the family of orbits with the same angular momentum. Consider the 



effective potentials of figures 8(a), and 9(a), which show V e fr's ranging from large L 
to low L, or equivalently successive Kff's from the weak-field to the strong- field regimes. 
Clearly, the stable circular orbit always has the lowest energy for a given L. If, as we will 
show, A(p r = 27r(g + 1) increases monotonically with energy, it must also be the case that 
q = w + v/z increases monotonically. It follows immediately that eqn. (|T2j) is true. We 
verify these claims for the Schwarzschild spacetime in this section. 

Spinning black hole spacetimes do not admit a comparable one-dimensional effective 
potential description. However, all the qualitative features of the Schwarzschild analysis 
survive and, in fact, all the interesting dynamical effects become more pronounced. We 
outline the dynamical results for the spinning black hole spacetime in §111 Bl 

We begin in the weak-field regime at large radii (large L) and will watch the dynamics 
evolve as we move into the strong- field regime at close separations (low L). Consider, for 
instance, solar system values, so that L is in the range of Mercury's angular momentum 
L ~ 2 x 10 4 . Such a potential at large r is shown in figure [71 There is effectively a 
large centrifugal barrier at low r. The energy of the bound orbits ranges from E = 1 for 
the eccentricity e = 1 trajectory down to the stable circular orbit e = 0. (Following the 
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standard convention in the literature we define eccentricity as 

r„ - r„ 



(15) 



r a + r p 

where r a is the apastron and r p is the perihelion.) 

All of the bound orbits with this L are well described as minute precessions of an ellipse. 
In our taxonmy, these orbits correspond to (z, 0, 1) with a very high number of zooms z 
(q = w + v/z << 1). Mercury, for instance, precesses roughly 43"/century. With an orbital 
period of 88 days, this comes to roughly 0.1"/orbit, which is ~ 360°/(3600 2 ). Mercury, 
therefore, is close to a periodic orbit with z = 3600 2 = 1.296 x 10 7 leaves and no whirls, 
which always advances to the next available apastron - that is, a (1.296 x 10 7 , 0, 1). Even the 
perfectly periodic (1.296 x 10 7 , 0, 1) will give the appearance of a simple precessing ellipse. 

In figure [HI we move into the intermmediate-field regime and show the potentials for 
angular momenta L = 4.1 and L = 4.5, with the centrifugal barrier better resolved to show 
the unstable E > 1 circular orbit. 



The key conclusions for the black hole dynamics comes from the simple figures 8(b) and 



8(c) Some salient features turn out to be universal for all L. Most important, figures 8(b) 



and 8(c) show that, for a fixed L, q = w + vj z increases monotonically with both E and 
e and terminates at some finite value. As per the comments at the beginning of §111} this 
justifies eqn. ([121 . 

There are other interesting dynamical features evident from these figures alone. For 
instance, we see that no orbits in this intermediate L range whirl since q = w + v/z<l 
and therefore w = over the entire range of eccentricities. Incidentally, the same thing 
is true for the planetary L (whose q versus e plot we omitted because it was graphically 
indistinguishable from the q = axis). This explains why we do not see zoom- whirl behavior 
in the solar system, or for that matter, even for Schwarzschild orbits whose periastra are as 
close to the black hole as r « 8. 

For L = 4.5, the allowed closed orbits correspond to rational numbers bounded by the 
range 

1 3 

, (16) 

where we have approximated the limits by the nearest low v, z option (w = 0). It follows that 
all orbits for L = 4.5 will look like precessions of (4, 0, 1) not quite going over to precessions 
of (3, 0, 1) at the upper limit as they march through the rationals with energy increasing to 
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(a)Thc effective radial potential V e ^[r) as a function of r on a log 
scale for L just above Ljbco- The horizontal line at V e Q =1/2 
corresponds to particle energy E = 1. The y-axis ranges from 
0.5 + 9 x 1CT 9 down to 0.5 - 3 x 10~ 9 . 



Schwarzschild 
(o = 0) 



Schwarzschild 
(o = 0) 




(b)w + v/z versus energy. (c)w + v/z versus eccentricity. 

FIG. 8: Schwarzschild, values of L just above Ljbco- I n terms of accumulated angle, (A(^ r /27r) 
1 = w + v/z. 
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E = 1. This is dramatically different from Mercury- type motion. These are not precessing 
ellipses but are precessions of a four-leaf clover. 

As the angular momentum drops even further to L — 4.1, the allowed closed orbits 
correspond to rational numbers bounded by the range 



- < q < g , (17) 



where we have approximated the limits by the nearest low v, z option. Orbits will pass 
through 3/8, 2/5, 3/7, 1/2, 4/7, to 3/5 and all the rationals in between as the energy increases 
from that of the stable circular to E = 1. 

The zoo becomes more crowded as the angular momentum drops to the point that the 



energy of the unstable circular orbit drops below 1. The top darkened curve in figure 9(a^ 



is the potential for L = Ljbco — 4, the angular momentum of the innermost bound circular 
orbit (IBCO). The bottom darkened curve is the potential for L = L IS co — VT2, the 
angular momentum of the innermost stable circular orbit (ISCO). For every Ljsco < L < 



Ljbcoi several of which figure 9(a) shows, there is an energetically bound (E < 1) unstable 
circular orbit. Accompanying each bound unstable circular orbit is an eccentric orbit called 
a homoclinic orbit. These orbits are so central to the discussion of strong- field dynamics 
that we now briefly outline how they fit into the (z, w, v) taxonomy before continuing with 
the rest of the dynamical discussion. 

A homoclinic orbit asymptotically recedes from a periodic orbit (in this case, an unstable 
circular orbit) in the infinite past, zooms out to some finite apastron, and then asymptotically 
approaches the same periodic orbit in the infinite future. Though eccentric, each homoclinic 
orbit has the same energy E u and angular momentum L u as the unstable circular orbit to 
which it asymptotes. When it exists, a homoclinic orbit forms the boundary (for a given L) 
between orbits that plunge over the top of the potential (E > E u ) and orbits that do not 
(E < E u ). 

To see the role of homoclinic orbits in our taxonomy, consider the behavior of a homoclinic 
orbit that starts at apastron, such as the one in figure [TUJ As the particle climbs toward the 
maximum of the effective potential, it asymptotes toward the circular orbit located there and 
executes an infinite number of whirls. Because it never returns to apastron, the homoclinic 
orbit is not strictly periodic. Nevertheless, it is appropriate to think of this single leaf orbit 
as the w — > oo limit of the progression of the (l,w,0) orbits. We thus label each homoclinic 
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(a)The effective radial potential V e &(r) bounded on the top by 
Libco = 4 and on the bottom by Ljsco = vl2. 



Schwarzschild 
(a = 0) 



Schwarzschild 
(o = 0) 




0.95 0.96 0.97 0.98 0.99 



(b)u> + v I z versus energy. 



(c)w + v 1 1 z versus eccentricity. 

FIG. 9: Schwarzschild for Ljsco < L < Libco- I n terms of accumulated angle, (A(p r /2w) — 1 = 
w + v/z. 

orbit as (1, oo, 0), so that each has an associated rational q = oo. 

It is perhaps surprising that the q = oo orbits are bound and have eccentricities e < 1. 
More specifically, homoclinic orbits with L — > Libco have eccentricities e — > 1, while those 
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FIG. 10: A homoclinic orbit (1, oo,0) for Ljsco < L < Ljbco- The orbital parameters are 
a = 0, L = 3.636619, and E = 0.958373. 

with L — > Ljsco have e — > 0. In fact, the ISCO is just the zero eccentricity homoclinic orbit. 
Properties of homoclinic orbits beyond what is needed for the current discussion, including 
analytic expressions for their trajectories, can be found in Refs. 13j, 121 ]. 

Since the homoclinic orbit is the maximum energy bound orbit (for a given L, it defines 
the upper range of allowed periodic orbits, a la eqn. (112j) . Therefore q ma x = oo for all 
Li sco < L < Ljbco- As always, the lower bound is greater than zero because q terminates 
at a finite value. 



To see the consequences of this fact, consider the L = 3.9 line in figure 9(b) , which shows 
a q c of approximately 0.475 ~ 19/40. By eqn. ( |T2l) . we have the bound 
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— < q < oo 
40 ~ H ~ 



(18) 



for L = 3.9. We can approximate the lower limit better by going to higher v, z, but these 
values are sufficient to make a point: there are no Mercury type motions. Rather, small 
perturbations of the circular orbit can be approximated by a periodic orbit with 40 leaves 
that jumps 19 apastra away every radial cycle. Such an orbit will look very much like a 
precession of a two- leaf orbit. Again, the one-leaf, zero- whirl orbit is forbidden. The (1,0,0) 



corresponds to a q = and is not reached by any L = 3.9 orbit as figure 9(b) shows 



Finally, figure [TT] shows a periodic table of orbits (still for L = 3.9). For illustration, we 
only inlcude orbits with z up to 3, although the spectrum allows z — > oo orbits. The first 
column corresponds to w — 0, the middle column to w — 1 and the last column to w — 2. 
The series could be continued through to w — > oo. The sequence should be read from top 
to bottom and then from left to right to indicate increasing energy. Notice that the first 
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FIG. 11: All z = 1,2,3 orbits with w = for the first column, w = 1 for the middle column, and 
w = 2 for the last column. All orbits lie on the L = 3.9 line of figure [9J Norbits increase in energy 
from top to bottom and left to right. The first radial cycle is emphasized in bold for each orbit. 
Notice the first and second entry in the first column are blank, indicating the inaccessibility of the 
(1,0,0) and (3,0,1) orbits. 
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FIG. 12: A series of w = 1 orbits for the L = 3.9 line of figure [9l Orbits increase in energy from top 
to bottom and left to right. All z = 1,2,3,4 orbits are shown. Also shown are randomly selected 
high z orbits. Notice that the high z orbits look like precessions of the energetically closest low z 
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entry in the sequence of periodic orbits is blank to indicate the complete absence of the 
(1, 0, 0) orbit as discussed above. The second entry in the sequence is also blank indicating 
the complete absence of the (3, 0, 1) orbit - as figure [9] shows, q never drops to 1/3. 

Figure fT2l which was previewed in the introduction, shows the progression through orbits 
with higher numbers of zooms. All orbits are drawn from the w — 1 band - the middle 
column of figure [TTJ As before they are arranged in order of increasing energy from top to 
bottom and left to right. All orbits with z — 1, 2, 3, 4 are drawn. Between each of these low 
leaf orbits, randomly selected high zoom orbits are shown as well. The high zoom orbits 
look like precessions of the low zoom orbits. 

Notice that for Schwarzschild the vast majority of low leaf orbits are stacked at high 
eccentricity. Orbits around spinning black holes, by contrast, show a wider spread of low z 
orbits at low eccentricities as we now demonstrate. 



B. The Kerr zoo 

Kerr orbits do not admit a simple one-dimensional effective potential description. 6 De- 
spite this complication, the generic Kerr system (appendix §|A~)) has all the key features we 
need for our taxonomy. For any L, there is always one stable circular orbit. At some critical 
L a first bound unstable circular orbit will appear bringing with it an associated homoclinic 
orbit. The unstable circular orbit always has energy higher than that of the stable circular 
orbit. For a given L IS co < L < L IBC o, the energy continues to specify uniquely orbits 
ranging from the stable circular to the homoclinc. Therefore, the same graphical analysis 
that we used to determine the population of the Schwarzschild zoo can be used to determine 
the population of the Kerr zoo. 

Because the weak-field regimes look essentially the same for all a, we consider only values 
of Lisco < L < Ljbco- As demonstrated for a Kerr black hole with spin a = 0.995 in 



figure 13(a), the rational numbers corresponding to closed orbits increase monotonically 



with energy. The rationals diverge at the homoclinic orbit for each L since the number of 



6 There is an unappetizing effective potential description for equatorial orbits that we do not want to rely 
upon here. 
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(&)w + v/z versus energy. 
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(b)w + v/z versus eccentricity. 

FIG. 13: Prograde orbits around a spinning black hole with a = 0.995. The lines indicate increasing 
L from left to right through the values L = 1.57, 1.61, 1.82,2,2.1 



whirls diverges. The rationals are bounded from below by q c at the stable circular orbit: 



q c <q<oo 



(19) 



Figure 13(b) shows the increase of the rationals from the stable circular orbit to the 



homoclinic orbit is shown versus eccentricity for a = 0.995 and a range of L's. 

In contrast to the non-spinning case of figure M, all eccentric orbits (with L < Ljbco) 
have at least one whirl for such high spin. Generally, a ^ orbits exhibit more whirls at low 
and intermediate eccentricities than their a = counterparts. As eccentricity (and energy) 
increase for a given L, the number of whirls diverges as a homoclinic orbit is approached. 

Figure [14] displays a periodic table, energy increasing from top to bottom through the 
w = 3 and w = 4 bands for L = 1.82. Each entry has a periodic orbit. All z = 1,2,3,4 
closed orbits are shown. The series shows the passage through these low leaf, moderate 
eccentricity orbits as energy - and therefore q increases. 

Finally, figure [15] displays a periodic table with both low and high z orbits, energy in- 
creasing from top to bottom for L = 2 and a = 0.995. The sequence takes you through low 
leaf orbits to high leaf orbits that look very much like precessions of lower leaf periodics. 
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FIG. 14: A series of w = 3 and w = 4 orbits for the L = 1.82 line of figure 1 131 Orbits increase in 
energy from top to bottom and left to right. All z = 1, 2, 3, 4 orbits are shown. 
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FIG. 15: A series of w = 3 orbits for the L = 2 line of figure PT31 Orbits increase in energy from top 
to bottom and left to right. All z = 1,2,3,4 orbits are shown. Also shown are randomly selected 
high z orbits. Notice that the high z orbits look like precessions of the energetically closest low z 
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The conclusion is that all eccentric Kerr orbits show zoom-whirl behavior of some kind 
for this a and L < Ljbco- m the strong-field regime, once L < Libco, the pattern of 
zoom-whirls at the lower energy bound looks like a precession of a low-leaf clover and at the 
upper energy bound marches toward homoclinic - a single leaf with an infinite number of 
whirls. 

RECAP OF DYNAMICAL RESULTS 

This concludes the primary results of the paper. Before turning to a discussion of astro- 
physical applications of the taxonomy, we briefly recap the results of this section: 

• For a given a and L, all periodic orbits corresponding to rationals in the range 

q c <q< qmax (20) 

define the skeleton of the entire orbital dynamics, where q c is the limiting rational associated 
with the stable circular orbit and q ma x is the approximate rational of the maximum energy 
bound orbit. 

• In the Newtonian limit, the upper and lower bounds both approach zero; Keplerian 
orbits are ellipses. 

• There is no zoom-whirl behavior in the weak-field regime. 

• At the ISCO, the upper and lower q bounds both approach oo because the ISCO is at 
once both a circular and a homoclinic orbit and all homoclinic orbits have q = oo. 

• In the very strong-field regime, which we take here to correspond to (L < L IBC o), q 
has no upper bound: 

q c < q < oo (21) 

• In the strong-field regime, the simple precessing ellipse familiar from planetary orbits 
is forbidden. 

• All aperiodic eccentric orbits are precessions of low-leaf periodics. 

• As the ISCO is approached, all orbits whirl as well as zoom for any a. 

IV. APPLICATIONS AND FUTURE WORK 

Our approach - representing the entire black hole dynamics in terms of a periodic skeleton 
- can be extended to non-equatorial orbits and then applied to astrophysical problems such 
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as gravitational wave astronomy and a dynamical analysis of spinning black hole pairs. We 
briefly describe each of these applications in turn. 

Non-equatorial Orbits: Equatorial orbits yield especially well to a visual rendering of 
the properties of periodic orbits and their relationship to aperiodic ones. We focused here 
exclusively on equatorial orbits for precisely those pedagogic reasons. Already underway, 
our next order of business is to extend the periodic orbit analysis to non-equatorial Kerr 
geodesies. Generic non-equatorial trajectories can look complex particularly for high eccen- 
tricities and incliniations in the strong-field [lij. Still, periodic orbits must exist. Just as 
in the equatorial case periodic orbits will correspond to rational numbers. Identifying the 
ordering of those rationals with respect to those orbital parameters should offer an elegant 
and illuminating description of the dynamics. 

Framing the non-equatorial periodic skeleton is not a solely academic enterprise. The 
modeling of extreme mass ratio inspirals (EMRIs) of compact objects into rotating super- 
massive black holes relies heavily on non-equatorial Kerr dynamics. 

Gravitational Wave Astronomy: Compact objects orbiting supermassive black holes 
in galactic centers are expected to be plentiful sources for gravitational waves in the LISA 



bandwidth 
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161 ] . Tens of thousands of such compact objects populate our own galactic 



center. Through gravitational scattering, a stellar mass compact object can be thrown into 



a highly eccentrc orbit in the strong- field regime [26j. We have shown that all eccentric 
orbits show zoom behavior of some sort - that is to say, a pattern of discernible leaves - 
and large ranges of L and E generate whirl behavior (figures 191 and fT3i) . In particular, in the 
strong-field regime, we should not expect eccentric orbits to be characterized by precessing 
ellipses, but rather by the precession of multi-leaf orbits. These EMRIs will therefore show 
the zoom-whirl behavior taxonomized in this paper. 

So far we have only described the conservative dynamical system. Astrophysical objects 
orbiting near a massive black hole will radiate gravitational waves and veer off geodesies 
under the effects of dissipation. A study of adiabatic inspiral through periodic orbits might 
offer computational advantages. For instance, periodic orbits should have more rapidly 
converging Fourier series, than aperiodic ones. Additionally, the evolution of the conserved 
quantities E, L, and Q may not be independent, thereby cutting down dramatically on the 
computational expenses. 

Closed orbits also show promise as a benefit to signal extraction from the gravitational 
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wave observatories in a time-frequencey analysis. Furthermore, the templates of aperiodic 
orbits, where required [2^], might be suitably replaced by periodic tables. 

Spinning Pairs: The Kerr orbits are only integrable in the case of non-spinning test 
particles. If the companion spins, a constant of integration is lost (namely Q) and the 
dynamics is no longer integrable 28|, [29] . While many of the other techniques for studying 
the dynamics may not generalize to spinning pairs, periodic orbits still define the skeleton of 
the phase space. Even if no clean taxonomy for those orbits exists, a search for and analysis 
of periodic orbits should reveal general properties of the dynamics. 
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urther, it is already known that spinning pairs can exhibit chaos 28 
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35( . Treating spin of the companion as a perturbation to Kerr motion, the progression of 



closed oribts under that perturbation could track the transition to chaos. As the companion 
spin increases, the periodic orbits proliferate and the homoclinic orbit is replaced with a 
homoclinic tangle - a fractal set of unstable periodic orbits Q|- A fractal in phase space 
of such rapidly spinning systems was revealed in Ref. 29J for comparable mass black hole 
binaries in the Post-Newtonian expansion by numerically scattering black holes. An analysis 
centered on periodic orbits should reveal the fractal set more readily, thereby removing the 
need to search and survey phase space blindly for chaotic behavior. 



V. CONCLUSION 



Periodic orbits structure the entire equatorial dynamics around black holes. In the spirit 
of Poincare's approach, we have defined a taxonomy of all periodic orbits based on zooms, 
whirls, and vertices and forged an explicit connection with the set of rational numbers. 
Under this scheme, all generic aperiodic orbits can always be approximated as near one of 
the periodics. 

Our method reveals that no eccentric orbits will trace out planetary type elliptical pre- 
cessions in the strong-field regime. Instead, all eccentric orbits are precessions of multi-leaf 
orbits with whirls multiplying as the black hole spins. The implication for gravitational 
wave astronomy is intriguing: all eccentric orbits will transition through these nearly peri- 
odic orbits as they inspiral to the plunge. 
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APPENDIX A: THE KERR EQUATIONS 

We now present the explicit equations integrated to generate the results for §111 Al and 
IIII Bl We will also derive the association between periodic orbits and rational numbers 
from the dynamical systems perspective. Since the natural language for that discussion is 
Hamiltonian mechanics, we begin with the Hamiltonian formulation of Kerr geodesic motion. 



1. The Kerr Metric and Geodesic equations 

In Boyer-Lindquist coordinates and geometrized units (G = c = 1), the Kerr metric is 

,2 2,2 ( 2Mr\ , 2 AMar sin 2 6 , , S , a „ in2 
ds 2 = -/j, 2 dT 2 = - ( 1 — J dt 2 dtdcp + —dr 2 + Ed6 2 

. 2n ( 2 2 2Ma 2 rsin 2 6\ , 2 , Ai x 
+ sin 2 6 [r 2 + a 2 H dy? 2 . (Al) 



Here 

S = r 2 + a 2 cos 2 

A = r 2 - 2Mr + a 2 , (A2) 

and M, a denote the central black hole mass and spin angular momentum per unit mass, 
respectively. Besides the test particle mass /x, three additional quantities are conserved along 
Kerr geodesies: the energy E and ^-component of angular momentum L as measured by 
observers at infinity, and the Carter constant Q. Following a useful and now common con- 
vention, we hereafter set both M — 1 and \i — 1. This choice leaves the overall appearance 
of the metric and other Kerr expressions otherwise unchanged provided all Boyer-Lindquist 
coordinates, the proper time r, the spin parameter a and all conserved quantities are now 
interpreted as dimensionless quantities. 

Since they possess as many constants of motion (/i, E, L, Q) as degrees of freedom 
(t,r,9,<f), the usually second order geodesic equations can be integrated to yield a set 
of 4 first order equations of motion for the coordinates 4J, 



Sr = 


±^ 


IR 


£fl = 


±-\ 


/e 


Htp = 




aE 



L \ a 
^9 + A P 
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r 2 + a 2 



Hi = -a (aE sin 2 6 - L) + ^ P , (A3) 

where an overdot denotes differentiation with respect to the particle's proper time r and 

f T 2 



9(0) = Q-cos z 9<a 2 (l-E , 

[ sin 6 1 

P(r) = £(r 2 + a 2 ) - aL 

R{r) = P 2 - A {r 2 + (L — a£) 2 + Q} . (A4) 

We note that while these equations are concise and appealing in some ways, during 
numerical integration they tend to accumulate error at the turning points due to the explicit 
square roots in the r and 9 equations, not to mention the nuisance of having to change the 
signs of the r and 9 velocities by hand at every turning point. While other authors circumvent 
this problem with a rather involved reparametrization of the equations in ( 1A3I) . we will see 
that in a Hamiltonian formulation the numerical difficulties are avoided naturally. 



2. Hamiltonian Formulation of Kerr Geodesic Motion 

We begin with a Lagrangian for a free particle in the Kerr spacetime: 

£ = g&tf ?V • ( A5 ) 

Because the q a are dimensionless coordinates (fi = 1), and because Qa^q 13 = — 1 for timelike 
trajectories, eqn. flAlj) implies that C is identically equal to —1/2 along any trajectory. The 
particle's 4-momentum 

P a = q a 

is also dimensionless and is identical to the 4- velocity. Defining the canonical momentum p a 
for the system in the standard way, we see that it coincides with the 4-momentum one-form, 

* s w (A6) 

= 9*p(f = 9af3P^ ■ (A7) 
7 



Explicitly, the components of the momentum are 



2r \ ■ 2ar sin 2 9 . . , . 

p t = - (1-— Jt cp (A8) 



7 "Momentum" will always refer to the canonical momentum, i.e. the momentum one-form on the spacetime, 
unless otherwise noted. 
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Pr = -£r (A9) 

p e = T.9 (A10) 

■ 2 n ( i i 2a 2 r sin 2 9\ . . . 

p v = sin^ ( r 2 + a 2 + J ip (All) 

2ar sin 2 6* . 

■t . (A12) 



E 

We finally define the Hamiltonian in the standard way, 

H = p^-C= l -g^p a p p . (A13) 

Note that the Hamiltonian is numerically identical to our Lagrangian (—1/2) along any orbit 
because each quantity is just half the contraction of the 4-momentum. This is not surprising 
since, for geodesic motion, the Lagrangian and Hamiltonian contain only kinetic terms. 

To get an explicit expression for the Hamiltonian in terms of our dimensionless coordinates 
and momenta, we could compute the inverse metric g a @ and grind through eqn. ( 1A13I) by 
brute force. However, a simple observation yields a nice expression for H(q,p) with much 
less effort. Using eqns. flA9l) - f)A10j) in eqn. flA5j) . we know the Hamiltonian must have the 
form 

H(q,p) = ^p 2 r + ^p 2 e + f(r,6,p u p^ , (A14) 

where / is a function yet to be determined that depends only on p t and Plf and, via the g a ", 
on r and 9. Eliminating the p r and pg above with the use of (1A8I) and (1A3I) and noting that 
H must identically equal —1/2, we see that 

H(q, P) = ^ + ^ + f(r, 9, p t , Plp ) = -1/2 . (A15) 

This fixes / and allows us finally to write H as 

tt, \ A 2 1 2 R+AQ 1 , Aie , 
*(*P) = ffirf + ffirf--W--2 ' (A16) 

where R and G are the functions listed in eqns. ( 1A4I) with every E or L therein interpreted 
as a — p t or P(p , respectively, and with every Q treated as a constant. 
With the Hamiltonian written as in (1A16[) . Hamilton's equations 

• 9H . dH , Aify . 
Qi = jr- , Pi= -jr- (A17) 

dpi dqi 
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for the test particle motion become, explicitly. 

A 



s Pr (A18) 
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2AE dL 
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1 Y 2 fR + AQ 

■ 2 + 



2S/ ^ y V 2AS 



(A19) 



(R+ A6) 



(A20) 



where the superscripts ' and ^ denote differentiation with respect to r and 6, respectively. 
Note that since none of the equations contain square roots any more, they constitute a 
smoothly differentiable system of ODEs, even at turning points, and they can be integrated 
directly without resorting to an intricate change of variable. Since the p t and p' v equations 
vanish, the only added cost of retaining analytical transparency in the equations is that we 
must integrate 6 equations in (1A18I) instead of 4 in (1A3I) . These are the quations we integrate 
numerically to identify the periodic orbits. 

We isolate the periodic orbits by inputing (z,w,v) for a given a and L to locate the E, 
apastron r a and perihelion r p of the corresponding periodic orbit. Numerically the periodics 
are extracted by explicitly integrating 

f t[ra) , f ra dip , / \ \ 

A<p r = 2 -f-dt = 2 -j-dr (A21) 



Jt{r p ) dt J rp dr 

to find those orbits for which the integral is a rational multiple of 2ir. 

As a final check of our numerical results, we run an independent code that computes the 
orbits directly from the geodesic equation and find the two independent computations in 
complete accord. 
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APPENDIX B: THE ACTION- ANGLE PICTURE 

In addition to deriving numerically supple equations of motion, we can use Hamiltonian 
dynamics to relate the periodic orbits to rationally related canonical frequencies. We now 
relate those frequencies to the taxonomy we have developed. 

It is a classic result of Hamiltonian dynamics that motion in a system with N coordinates 
and N conserved momenta will be confined to N- dimensional tori. In this case we do have 4 
coordinates (t,r,9,ip) and 4 conserved quantitites (//, E, L,Q). However, this classic result 
only applies to orbits bound in phase space. In our relativistic system, the t direction throws 
a wrench in the works since every orbit increases without bound in the timelike direction - 
worldlines, in both phase space and configuration space, are not compact. To circumvent 
this problem, we work in a formally reduced phase space. 8 Any canonical coordinate can be 
chosen as the time parameter and its canonical momentum will become the Hamiltonian for 



the reduced phase space 22j. In our case, we choose to parameterize the orbits by t and use 
Pt = E as the Hamiltonian of the reduced phase space. 

In this paper we are only dealing with equatorial motion, for which motion lies in the 
4D hypersurface defined by 9 = 7i/2,pg = and on which Q = 0. On that hypersurface 
r,ip,p r and p v still form a canonical set, so the hypersurface corresponding to equatorial 
orbits constitutes a bona fide 4D phase space. 

On this space, all phase space trajectories corresponding to orbits that are bound and 
non-plunging in configuration space lie on compact hypersurfaces topologically equivalent 
to 2D tori. Thus, we can make the transformation to action-angle variables for those orbits, 
to which we now turn. 

A canonical transformation to action variables allows us to set each canonical momentum 
Jj equal to a function of the constants of motion. The reduced Hamiltonian can be rewritten 
as a function of the J only, 

E = E(J) , (Bl) 

so that it is cyclic in all the ipi, and the (ipi, Jj) form a canonical set whose evolution is still 
determined by Hamilton's equations 

d^h dE 
dt d^i 



A more detailed exposition of which we leave to the noncquatorial case. 
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dt 



dE 



:b2) 



(i = r,tp). Together, (1B1I) and (1B2|) imply that each Jj is a constant of the motion and that 
each ipi is linear in time, 

^=0^ + ^(0) , (B3) 



where each 



U)i(J) 



dE(J) 

dJi 



(B4) 



is a function only of the J, and thus also a constant of the motion. These constants are the 
orbital frequencies used as a basis in which to Fourier decompose functions of generic Kerr 
orbits, including the instantaneous adiabatic gravitational waveforms they emit 
We define the radial and azimuthal actions in the standard way, as 



23 



24]. 



J, 



— (h p r dr 
2~ 



(B5) 
(B6) 



where the curves C r and C v are projections of the orbit into the r — p r and ip — planes, 
respectively. Because p v is a constant we see that 



J, 



— Pip 



(B7) 



We are more interested in the frequencies u r and uj v . An explicit expression for the 
Hamiltonian in terms of the J's requires inverting the J r integral in (1B5H to solve for 
E(J r , L) = E(J r , J<p) and thus is not analytically accessible. Still, expressions for the Ui 
in terms of E, L rather than in terms of J r , J v can be computed as follows. 

Since a trajectory can alternately be specified by constants of the motion (E, L) or 
(Jr, Jip), there exists some transformation from the (E, L) set to the J r , J v set. Since 



J ip {E,L) = L 



;bs) 



the Jacobian of the transformation is 



dJr 


dJ r 




8J r 


OJr 


dE 


0L 




dE 


dL 


8J V 


dJtp 







1 


~dE 


~dL 





0J r 
'dE 



(B9) 
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which, since J r increases monotonically with E at a given L, never vanishes. There thus 
exists an inverse transformation E(J r , J v ), L(J r , J v ) from the actions back to the usual 
orbital constants. The product of those two transformations should be identity: 



dJr 8Jr \ I 8E_ 8E_ \ / dJr dJr \ / dE_ dE_ 

dE dL I I dJ r dJ v I _ I dE dL \ [ dJ r dJ v 



8L dL 



1 / \ 



dE dL 

which is equivalent to the two independent equations 

dJ r dE 




(BIO) 



dE dJ r 

dJ r dE dJ, 



1 (Bll) 

(B12) 



dE dJ v dL 

But in our reduced phase space, E is the reduced Hamiltonian, so each dE/dJi = u>i. The 
equations above then yield 

1 



dJ r /dE 
dJ r /dL 
dJ r /dE 



^ = --^T^ ■ (B13) 



Combining equations (1A8|) . (1A4I) and (1B5I) . we see that 



Or 



p r dr = 2 I dr^— (B14) 



A 



and, combining this with the equations of motion, that 



dL 2tt dL J r A 



1 [ ra 1 / dR 
-2 / dr 



2vr J rp 2 v / i?A V dL 

where the term proportional to the partial of the limits vanishes since the integrand vanishes 
at the limits. Continuing, 

dJ r 1 „ f ra , £ -dR/dL 
-2 / ar- 



dL 2tx J rp 2AE 
-2 / dr- 



2tx L r 

^-Acp r , (B15) 
271 
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where the last equality comes from equation (J7J). Additionally, 



»* J_ 2 » f r "dr VR 



dE 2n dE J r A 

1 f ra , 1 /&R 
-2 / rfr 



2tt J rp 2y/RA \dE 

1 r a , £ dR/dE 
-2 / ar- 



2tt J rp y/R 2AE 

i r a , i 

-2 / dr- 



2tt J r r 

= 7^T r , (B16) 
2ix 

where T r is the radial period of an orbit. 

The physical interpretation of the tu's is now clear. As expected, u r is 

2tt OR 2tt 
Ur = T r d^dE = ^ ' (B17) 

where T r as measured in coordinate time t. The rate at which azimuth accumulates averaged 
over the radial period is u v : 

Jrp ^ AVR 8L A(^ r 
Jr p ar AVRdE 

Because for a periodic orbit, the total orbital period is an integer multiple of T r (in fact, 
zT r ), we can equivalently say that uj^ is the rate of increase of tp averaged over the entire 
orbital period (this latter interpretation is what will generalize when we leave the equatorial 
plane) . 

Finally, to make the connection with the §11 B[ note that the frequency ratio 

^ = ^ = l + ^ = l + g . (B19) 

L0 r 2ix z 

Thus, in the action-angle picture, the periodic orbits are those with rationally related fun- 
damental frequencies. 

A typical orbit will fill out the torus because the canonical angular frequencies are not 
commensurate. In contrast, closed orbits do not. It is precisely this relation we have 
exploited to explicitly locate the periodic orbits in our taxonomy. 
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APPENDIX C: CIRCULAR ORBITS REVISITED 

As discussed briefly in §11 C\ circular orbits, while clearly periodic, are somewhat anoma- 
lous in the geometric picture and our (z, w, v) scheme. After all, our taxonomy hinges on 
quantities calculated per radial period. Since the radial period of a circular orbit is zero, it 
is not clear what rational number to associate with circular orbits. Nevertheless, it turns 
our that in a useful and enlightening sense, some stable circular orbits (a measure zero set 
of them, to be precise) are mappable to the rational numbers. Naively, we might expect 
them to correspond to the rational number 1, just like Keplerian circular orbits, but this 
turns out to be a mistake. 

The action-angle picture is significantly more informative here. Our difficulty stems from 
the fact that all circular orbits have a vanishing radial action J r . In a loose sense, the 
actions correspond to the two circumferences that characterize a two-dimensional torus. 
When one of those circumferences is zero, the 2-torus T 2 just becomes a 1-torus T 1 = S l : 
a circle. Orbits that are circular in configuration space thus live on surfaces in the phase 
space that also have the topology of circles. In that sense, circular orbits only really have 
one fundamental frequency, u^, and there is no rational frequency ratio to speak of. They 
are all nonetheless periodic, for any value of u^, because any curve that fills out S 1 with a 
constant velocity necessarily closes on itself and becomes periodic. 9 

Nonetheless some stable circular orbits can be mapped to the rational numbers. We 
can construct that map in two equivalent ways. First, for a given L, we can take the zero 
eccentricity limit of the values of ou r ,oUtp, and their ratio and define the corresponding values 
for the circular orbit with that same L to be those limits. This was the attitude taken in 
aiTCl 

Alternately, we can perform a linear stability analysis of the circular orbits. We find the 
frequencies of small oscillations of the r and <p motions of low eccentricity orbits around the r 

9 Note that for spherical orbits, i.e. constant r orbits not confined to the equatorial plane, which other 
authors usually refer to as "non-equatorial circular orbits" , this will no longer be the case. Non-equatorial 
orbits in general have 3 associated actions J ri Jg and J v and are confined to surfaces with the topology of 
T 3 . On spherical orbits, J r vanishes, but the remaining actions do not, leaving spherical orbits to occupy 
surfaces with the topology of T 2 . Much like generic equatorial orbits, then, spherical orbits will only be 
exactly periodic only when the ratio of their frequencies loq /uj v is rational. 
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FIG. 16: Left: versus the radii of stable circular orbits. Right: lo t versus the radii of stable 
circular orbits. Increasing from right to left the spin values are a = 0,0.5,0.8,0.95,0.995,0.998. 
All orbits are prograde. 

and ip motions of the reference circular orbits and examine the ratios of those frequencies 10 . 

Whether derived from a stability analysis or as the limiting values of low eccentricity 
orbits, figure [16] shows the resulting values of uo^ and uo r while figure [17] shows u) r for 
stable circular orbits as a function of their radial coordinate r for various values of a. Some 
features are worth noting. First, u v for a circular orbit is just the coordinate angular velocity 
and thus satisfies 



the well-known relativistic generalization of Kepler's third law [25|. Second, the uv's differ 
from the Keplerian in that the increase with decreasing r until they hit some maximum (in 
the Keplerian case, the u r continue to increase and diverge at r = 0), at which point they 
decrease, reaching zero at the ISCO. The ISCO, then, will have a diverging uJ v /u r ratio, an 
observation that was relevant in our discussion of homoclinic orbits. All homoclinic orbits 
have q = w + v / ' z = u v /uj r = oo, including the ISCO is (the eccentricity zero homoclinic 
orbit). 

Third, u^/ujr for circular orbits is a continuous and monotonic function of r, increasing 
from an asymptotic value of 1 in the r — > oo limit and diverging as r — > rjsco- The continuity 
of this ratio means that although a measure zero subset of the circular orbits corresponds 
to rational numbers or, equivalently, to (z,w,v) triplets, most do not. Nevertheless, even 



id 



Often called epicycle frequencies, or just epicycles for short 



10 




FIG. 17: The ratio o^/ctv, which equals 1 + w + v/z, as a function of the radii of stable circular 
orbits. Increasing from right to left the spin values are a = 0, 0.5, 0.8, 0.95, 0.995, 0.998. All orbits 
are prograde. 

the irrational circulars are arbitrarily close to some rational. We can thus characterize 
every circular orbit either exactly (the rational circulars) or approximately (the irrational 
circulars) by a (z,w,v) triplet. 

Finally, dissipative dynamical systems with multiple coupled frequencies, can sometimes 
attractors (or at least transient attractors). Those attractors typically correspond to ra- 
tionals q with low denominator fractional parts (low z). The 1 : 1 and 3 : 2 spin-orbit 
frequency ratios of the moon-earth and Mercury-sun systems, respectively, are solar system 
examples of where this sort of behavior takes place. The possibility that an orbit decaying 
adiabatically under radiation reaction might become trapped in a resonance en route to 
plunge is thus open. At the time of submission of this article, we are still investigating that 
possibility. 
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